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A method for producing an upper bound for all multipartite purification protocols is devised, 
based on knowing the optimal protocol for purifying bipartite states. When applied to a range of 
noise models, both local and correlated, the optimality of certain protocols can be demonstrated for 
a variety of graph and valence bond states. 
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I. INTRODUCTION 

Quantum states of many qubits are essential ingredi- 
ents in the functioning of quantum computers, and yet 
their properties, such as entanglement, are poorly un- 
derstood. Of particular interest are graph states, which 
provide a range of benefits for communication and cryp- 
tography [l| . In addition, they form a universal resource 
for quantum computation ^, ^ and enable computation 
in scenarios where the employed two-qubit gate is prob- 
abihstic However, any practical implementation 

will introduce noise to the system. The noise-induced er- 
rors need to be minimized, and corrected, before a prac- 
tical application of these states is considered. One way to 
achieve this is by purification, where many copies of the 
noisy state are combined to yield a single perfect copy. 

The concentration of entanglement in two-qubit sys- 
tems has yielded some important results in quantum in- 
formation, such as secure communication via privacy am- 
plification 0. Nevertheless, the detailed study of similar 
transformations in many-qubit systems has only recently 
commenced. The purification of two-qubit states was 
first examined in 0, Q, and the performance of these 
protocols is optimal if the operations are perfect. The ef- 
ficacy of the proposed protocols in the presence of noise 
was explored in 9] . Subsequently, the question of purify- 
ing multipartite states has arisen [lO| . A variety of proto- 
cols have been discussed, starting from a subset of graph 
states [ni [13 and generalizing to arbitrary graph [l^ 
and stabilizer states [13] ■ These different protocols tend 
to trade between a large tolerance of noise U^,^ and the 
scaling of purification rate To date, little has been 
said on the subject of what is optimal, although some 
(non-tight) bounds have previously been found, such as 
in the case of independent local Z-noise [l^ , or for GHZ 

states 

In this paper, we prove the optimality of certain purifi- 
cation protocols for a variety of noise models by consider- 
ing general upper bounds. These derivations are based on 
a central theorem that analyzes the purification of mul- 
tipartite states in terms of the purification of bipartite 
states, and hence allows for a direct extension of previ- 
ous optimality proofs. When restricted to Z-noise, the 
application of the theorem becomes straightforward for a 
certain class of graph states called locally reconstructible 
states. These states allow for the direct application of 



the optimal bipartite purification protocol for each link 
of the graph, thus extending it to the multipartite case. 
Subsequently, a wide range of upper bounds in the toler- 
ated error rates are derived for a variety of error models 
and states, while optimality is numerically demonstrated 
in certain cases. 

The paper is organized as follows. After an initial 
introduction to graph states (Sec. [lT| and purification 
protocols (Sec. IIIip . we introduce our main theorem in 
Sec. nil CI which proves that if purification of a bipartite 
state is impossible, so is the purification of related multi- 
partite states. This result is applied to a variety of error 
models in Sec. IIVI including local Z-noise and depolar- 
izing noise. For the case of Z-noise we prove optimality 
of the two purification protocols under consideration for 
a sub-class of graph states. In the case of maximally 
depolarizing noise, we prove a universal bound which ap- 
plies to all graph states. Finally, in Sec. |Vl we derive an 
example of a valence bond state that can be optimally 
purified, proving that our method is not merely limited 
to graph states. Critically, the example that we produce 
has a finite entanglement length. The presented results 
expand and extend the work of [2(| . 



II. GRAPH STATES 

For the majority of this paper, we will be interested 
in the purification of graph states. These can be defined 
in two equivalent ways. With a particular graph G, we 
can associate a set of vertices, Vg, and edges, Eq, which 
connect pairs of vertices. The first way to define a pure 
graph state is as the ground state of the Haniiltonian 



H=j2 j^^^ n 
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where we have attached a qubit to each vertex and 
Xi and Zj are the familiar Pauli matrices applied to 
qubits i and j respectively. The individual terms Ki — 
Xi J^|,- j^^^Ea commute with each other, [AT^, Kj\ — 0, 
and thus stabilize the graph state. These are scaled by ar- 
bitrary coupling strengths Ji , which we will subsequently 
take to be equal. We note that [Ki, Zj] — if i ^ j, and 
{Ki, Zi} — 0, which means that the excited states of the 
Hamiltonian are described by local Z-rotations, and that 
these local rotations constitute a complete, orthogonal 
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FIG. 1: (a) A pentagon is locally reconstructible, but not 
two-colorable. This graph is locally equivalent to the five- 
qubit error correcting code, (b) The pentagon can be formed 
from two-qubit pairs by applying the projector P = |0) {00| + 
|1) {11| to each pair of qubits that needs to be combined, (c) 
The graph state which is locally equivalent to the codewords 
of the Steane-[[7,1,3]] error correcting code, (d) A shape which 
is two-colorable, but not locally reconstructible. 



basis over the Hilbert space associated with the graph. 
An equivalent way to define a graph state is in a more 
operational sense, where we prepare each qubit in the 
state |-|-) = (|0) -I- |l))/\/2, and apply controUed-phase 
gates along each edge of the graph. 

For what follows, the action of measurements on 
a graph state are of particular interest [llj. A Z- 
measurement removes a vertex from the graph, along 
with its edges. If we wish to form a graph state from 
two-qubit states corresponding to the edges of the graph, 
it is simplest to revert to the matrix product state for- 
malism [22], where we apply measurements on all of the 
qubits that need to be combined to a single qubit. In 
the case of a linear cluster state, this simply corresponds 
to performing a controUed-phase gate between the two 
pairs, measuring one qubit in the X-basis, and perform- 
ing a Hadamard rotation on the other. Although this 
simple method does not generalize to other graphs, the 
examples which we choose to give will be for linear clus- 
ter states, and hence this description is valid. For other 
graph states, the operation that we apply must project 
N qubits on a local site to 1 qubit, and is represented as 



P=|0) (0| 



For example, the required operations to create the pen- 
tagon of Fig. [TJa) are shown in Fig. [TJb). 



A. Locally Reconstructible States 

For a certain error model that we shall address (lo- 
cal .Z-noise), we will be particularly interested in the re- 
striction to a class of graph states which we call Locally 
Reconstructible (LR), defined as follows: 



Definition 1. Locally Reconstructible graph states are 
connected graph states for which there exists a non-trivial 
partitioning of the qubits into two parties such that nei- 
ther party has more than one edge from each qubit cross- 
ing the partition. 

In order to prove that this class of states is non-trivial, 
we examine some of its properties in App. |^ Impor- 
tantly, this class includes all cluster states (d-dimensional 
cubes), GHZ states (one vertex with edges to all others), 
and graphs which are locally equivalent to the codewords 
of error-correcting codes such as the Shor-[[9,l,3]] code, 
the 5-qubit code (Fig. [T](a)) and the Steane-[[7,1,3]] code 
(Fig. [IKc)). Indeed, in both of these figures, we provide 
a partition that demonstrates the LR character of the 
corresponding graph. 

We are not aware of this classification of graph states 
coinciding with any previous definition. For example, in 
Fig. [1] we provide two examples that show firstly that LR 
states are not necessarily two-colorable [ll] (Fig. [Ua)) 
and, secondly, that not all two-colorable states are LR 
(Fig. md)). In this case, for every possible partition, 
there is always a qubit that has at least two edges crossing 
the partition. One such partition is depicted in Fig.[Tl[d). 
Nevertheless, one can show that for all graphs of up to 
seven qubits the corresponding graph states are lo cally 
equivalent to LR states. Indeed, all the graphs in [2l| . 
which are used to categorize the local unitarily equiva- 
lent graphs of up to seven qubits, are LR. That is, one 
can reversibly transform pure graph states of seven or 
fewer qubits to LR states using local operations \/Ki. For 
larger systems, there exist examples that are not locally 
equivalent to an LR graph. Specifically, the icosahedral 
graph (12 vertices, degree 5, Fig. ^ and its local equiv- 
alences form a collection of 54 graphs, none of which are 
LR. Since the operations ^/TCi are Clifford operators, it 
remains a possibility that there are other local unitaries 
that could be applied which yield a different set of local 
equivalences [2ll. |23|. [2^ . 



III. PURIFICATION PROTOCOLS 

The aim of a purification protocol is to take many iden- 
tical copies of a noisy state pc and produce a single, pure, 
copy IV'g)- We consider that each qubit in the state pc is 
held by a different party (Alice, Bob. . . ), and that they 
hold the same qubit from every copy. The parties are 
restricted to applying Stochastic Local Operations and 
Classical Communication (SLOCC), which we initially 
assume to be perfect. These restrictions serve to illus- 
trate the entanglement properties of pa- There are also 
physically motivated systems where locality restrictions 
come into play, such as with quantum repeaters 9]. 

We assume that pc is diagonal in the graph state basis, 
i.e. 
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FIG. 2: The icosahedral graph, which is not equivalent to any 
LR states using the local operations \JKi. 



where j indicates which of the N qubits a Z-rotation 
is applied to. As has previously been explained else- 
where (see, for example, [13), any non-diagonal state 
can be made diagonal by SLOCC without changing the 
diagonal elements by probabilistically applying the sta- 
bilizers (with probabilities suitably chosen to negate the 
off-diagonal elements). The numbers < < 1 en- 
capsulate all the information about the errors. We are 
assuming that we know all of these values i.e. we know 
the noise model, although it is not strictly necessary to 
know the error probability, p. Hence, without loss of gen- 
erality, we assume that Ao is the largest element, since 
we wish to purify towards |'0g)- This is the fidelity of 
the initial state pa- If Ao were not the largest value, one 
could apply Z-rotations to make it the largest value, as 
these permute the diagonal elements. 

In this paper, we are interested in proving universally 
applicable bounds. We shall define the threshold fidelity 
/threshold ^s the valuc of Ao below which purification is 
impossible, regardless of the protocol employed, for a spe- 
cific noise model. All protocols have their own critical 
fidelity fcrit > /threshold bclow which that protocol does 
not work. Optimality of a protocol is proven when equal- 
ity holds. 



A. Genuine Multipartite Purification 

In the following subsection, we will propose a protocol 
that is based on bipartite purification and thus is easily 
analyzed. Indeed, for some types of noise and classes of 
states, one can prove its optimality. However, this proto- 
col will not perform well in general, and so when we are 
interested in other noise models, we are forced to resort 
to a multipartite_purification protocol, such as the one 
described in (l^ ■ We refer to this specific protocol as 
the Genuine Multipartite Purification Protocol (GMPP). 




FIG. 3: The Divide and Rebuild Purification Protocol. We 
start with many copies of the noisy graph state, (a) We form 
two-qubit nearest-neighbor states (noisy) by performing Z- 
measurements. (b) Two-qubit states are purified (if possible). 

(c) ControUed-phase gates are appfied between local qubits. 

(d) AU qubits except one from each party are measured in the 
X-basis, leaving the remaining qubits in the purified state. 

The GMPP is appfied to two-colorable graph states, 
which are graph states where every qubit can be labelled 
as either A or B such that all edges connect an A and a 
B. The protocol proceeds by application of arbitrarily or- 
dered sequences of two sub-protocols PI and P2. Since 
Z-rotations form a complete basis for the graph state, 
the state can be labelled by vectors fiA and /is, specify- 
ing which A or B qubits, respectively, have Z-rotations 
applied to them. The action of the two protocols is 

P"^ '■ ^HA.t^B ~ ^ ^. ^VA.jJ.B^IJ.ABl'A.tJ.B 
VA 

which then have to be renormalized. Both are realized by 
post-selecting on particular measurement results, which 
means that the rate of purification decreases exponen- 
tially with the number of qubits present. The application 
of the GMPP is challenging to analyze due to the arbi- 
trary choice of the sub-protocols PI and P2 at each step. 
In general, we resort to numerical exploration which, with 
finite computational resources, can never tell us precisely 
how close the GMPP comes to any upper-bound. An- 
other multipartite purification protocol has recently been 
proposed [l5j , which is easier to analyze for a range of er- 
rors, and achieves a superior purification rate. However, 
the critical fidelities at which it works are larger than for 
the GMPP, so little benefit can be derived from compar- 
ing them to the threshold fidelities which we calculate. 
Other recent work [l^ [3 has described protocols which 
are not limited to two-colorable states. 



B. Bipartite-based Purification 

The second purification protocol which we will analyze 
is certainly not new (see, for example, [i3j[iBi[lH|)j but 
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its simplicity enables the derivation of rigorous results. 
To implement the protocol (Fig. [3]), we initially measure 
the qubits of the graph state in the Z-basis until we are 
left with a single two-qubit state p2- Many copies of 
this state are then used to purify, if possible, a Bell pair 
\'4'2) ■ By performing Z-measurements on different sets of 
qubits, different Bell pairs are generated. Once we have 
a Bell pair for every edge in the graph, we can locally 
reconstruct the state, e.g. by applying controUed-phase 
gates and X-measurements. The conditions under which 
P2 — diag(Aon, Aqi, Aio, All) can be purified to \ip2) are 
well-known 0,[i], 

Aoo > 5, (2) 

and have been shown to be optimal using the positive 
partial transpose condition [H, US]- Hence, we only have 
to relate Aoo to the values of Aj of the original state, 
Pg- We refer to this protocol as the Divide and Rebuild 
Purification Protocol (DRPP). 

In [2(|, we considered the rate of purification for the 
DRPP, which is applicable whenever it is the optimal 
protocol. The rate of purification, R^, for the state jV') 
was described in terms of the rate of purification of a Bell 
state, i?2, taking the standard definition of rate, 

^ # of copies of 1-0) produced 
=?7^ of copies of p consumed 

This allowed us to bound the rate of purification by 

R2> Ri,> 

where iVgeo is a small geometric factor determined by the 
graph. This rate is a vast improvement over the GMPP, 
although other protocols have better rates at the cost of 
being less robust fT5l|. For d-dimensional cluster states, 
it was shown that the geometrical factor A^gco = 3c?^ is 
independent of the size of the graph. A similar argument 
can be applied to general graph states, and yields an 
upper bound 

TVgeo < min (^2{Dg - 1)Dg + 1, (^^^^ 

where Dg is the maximum degree of G i.e. no vertex 
has more than Dg edges. This does not coincide with 
the result for cluster states because we used knowledge 
of the geometry of cluster states to optimize our use of 
resources. 



C. Upper-Bound to the Purification of 
Multipartite States 

The intention of this paper is to make statements about 
when purification is impossible for all protocols. To pro- 
vide such a proof, we consider two parties Alice and Bob, 
each locally handling many qubits. The operations they 
perform are more general, but include, multipartite op- 
erations. 
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FIG. 4: If we assume the existence of a purification protocol 
for the multipartite state, then this implies that we can purify 
the two-qubit state, (a) Alice and Bob take the two-qubit 
state and reconstruct the noisy graph state, (b) This state is 
purified, (c) All extra qubits are measured out to return the 
original pair, now pure. 

Theorem 1. Consider the scenario where we wish to 
purify a two-qubit state p2. Provided that many copies of 
P2 can be converted by SLOCC into a noisy graph state 
Pg held by the two parties, then if purification of p2 is 
impossible, so is purification of pG- 

Proof. The veracity of this theorem is shown by contra- 
diction. We assume that purification of pG is possible 
by some protocol. From the condition of the theorem, 
we can start with p2 and locally convert it into pG by 
employing extra qubits. By assumption, there exists a 
protocol that can purify this state, leading to the pure 
state IV'g)- This can be converted into a maximally en- 
tangled two-qubit state by projecting out the additional 
qubits with local Z-measurements. Hence, p2 can always 
be purified if pG can be purified. If we know that p2 can- 
not be purified, we have a contradiction, and the initial 
assumption must be false. These steps are depicted in 

Fig.m 

□ 

A corollary is that if the state pG is described by a 
parameter p indicating the probability of an error occur- 
ring, this theorem gives an upper-bound on the value of 
p such that Pg can be purified. This could alternatively 
be viewed as a lower-bound on the required fidelity of the 
state Pg- We choose to refer to it as an upper-bound. 

Little is known about the conversion between bipartite 
(multi- level) mixed states pf " — *■ pG [11] , as required for 
the local reconstruction condition of Theorem [T] Thus, 
we have to examine different noise models on a case-by- 
case basis, which we do in Sec. lIVI For LR graphs, such as 
Fig. [3a), local reconstruction simply involves replacing 
each link across the Alice/Bob partition by a single copy 
of the two-qubit state. However, for non-LR graphs, such 
as Fig. [5]Jc) , this is not possible because two links need 
to be replaced which connect to a single qubit. While 
reconstruction is still possible in these cases, it generally 
means an increase in the local error probability, which 
becomes correlated in a different way to the errors in the 
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Q L> V 

FIG. 5: (a) Square graph shared between 4 parties, (b) Al- 
ice and Bob can locaUy reconstruct the square graph using 
two copies of p2, and applying controUed-phase gates between 
them, (c) The triangular graph is the simplest configuration 
for which the optimality proof fails. 

rest of the graph. This makes analysis more difficult, 
weakening the bounds that one can derive. 

Theorem [T] can be generalized by making two further 
observations. Firstly, it is not necessary to restrict to a 
bipartite state p2, since any state which can be locally 
converted into the state pc could be used. However, the 
only existing optimality conditions apply to two-qubit 
states. Secondly, the states that we use need not be graph 
states, it is just that the formalism of graph states guar- 
antees that we can convert \iPg) into \ip2)- In the general 
case, we should be able to perform measurements on lipc) 
which return a pure 2-qubit state with non-zero entangle- 
ment, and can subsequently be distilled to a maximally- 
entangled state. In Sec. |Vl we apply this method to the 
concrete example of a valence bond state. 

D. Comparison to the positive partial transpose 
condition 

Essential to the application of Theorem [1] is the knowl- 
edge of when a two-qubit state can be purified, which can 
only happen if there is non-zero distillable entanglement 
between the two parties. We can therefore interpret The- 
orem[I]as stating that purification of a multipartite state 
is impossible if there is a bipartite split for which there 
is no distillable entanglement. Following this interpre- 
tation, we can describe existing multipartite purification 
protocols as examples of bipartite purification protocols 
which use separable operations. In particular, we ob- 
serve that the operation of the GMPP is closely related to 
the protocol of [29] . The action of this high-dimensional 
bipartite purification protocol can be analyzed [30[, en- 
abling rigorous comparison between its performance and 
the upper-bounds calculated in this paper. 

With the aforementioned entanglement interpretation, 
it is simple to see that Theorem [T] is a constructive state- 
ment that can sometimes be applied to discover if a state 
can be written in the form 

i 

Consequently, we should expect similar results to those 
of where the positive partial transpose (PPT) con- 
dition is applied to bipartite groupings of multipartite 
states. However, in |19|, extra depolarization steps are 



applied which tend to remove entanglement from the sys- 
tem and, consequently, tight bounds are not expected. 

In general, how should our technique compare to the 
PPT condition, if one does not introduce additional de- 
polarizing steps? As expressed in terms of a recon- 
struction from two-qubit states, our upper bounds are 
strictly weaker than PPT. We can see that this is the 
case because starting from any two-qubit state which 
is separable (for which the PPT condition is necessary 
and sufhcient), and applying SLOCC results in a state 
which necessarily has PPT i.e. is separable, and purifica- 
tion is impossible. However, there also exist states with 
PPT which cannot be generated from two-qubit separa- 
ble states, which are known as bound-entangled states. 
Therefore our condition is strictly weaker than PPT. 

Nevertheless, our condition has two major benefits. 
Firstly, as already indicated, we need not be restricted 
to reconstructing from two-qubit states. In particular, if 
there exist bound entangled states with non-PPT (this 
still remains an open question) , then using these as a ba- 
sis for reconstruction yields a stronger bound on purifica- 
tion than the PPT condition can provide. Secondly, our 
technique is constructive, which eases its application in 
many scenarios, including the situation where the gates 
applied during the purification procedure are faulty (this 
will be explored in a later paper). 



IV. UPPER-BOUNDS FOR VARIOUS ERROR 
MODELS 

Given Theorem[l] it is interesting to apply the method 
to different types of noise, yielding bounds on when noisy 
states are not purifiable. It may not be possible to at- 
tain these bounds with purification protocols, and wc will 
be able to demonstrate that the DRPP docs not always 
achieve them. Numerical studies of the GMPP indicate 
a much tighter match in performance, although given the 
asymptotic approach to the bound, and the existence of 
strong local attractors, in most cases it is impossible to 
precisely verify whether they match. 



A. Local Z-Noise 

A straightforward application of our theorem comes 
when considering local Z-noise. While a very restrictive 
noise model, it has two physical motivations, namely that 
the thermal state of the Hamiltonian in Eqn. ([T]) is equiv- 
alent to the ground state with local Z-noise, and that Z- 
noise is a significant source of error in some experimental 
implementations, such as optical lattices. Moreover, as 
we will see, within our treatment this type of noise is the 
worst-case noise, giving the lowest probability threshold 
of all the local noise models considered here. 

We assume that a Z-error occurs on each qubit inde- 
pendently, with probability p. When we restrict to LR 
states, local reconstruction follows by replacing any links 
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across the bipartition with noisy two-qubit states. The 
structure of the class guarantees that the only subsequent 
operations that we need to perform are local controUed- 
phase gates between qubits. Since these gates commute 
with Z-errors, then starting with a two-qubit state, a 
many-qubit state can be built with the same error prob- 
ability. The two-qubit state cannot be purified if 

(1-P)'<i 

so this must hold for all LR states. 

Similarly, since Z-mcasurements commute with the Z- 
errors, we can show that the DRPP can purify the whole 
state provided 

Thus, the protocol is optimal, with a threshold probabil- 
ity of p = 1 - I/V2 ~ 0.29. Equivalently, 



/thr 



cshold 



= /crit = (l-p)" = 



1 



2^/2 



(3) 



This provides a useful benchmark to test other purifica- 
tion protocols, such as the GMPP (see Sec. IIV D|) . Note 
that the DRPP can purify all graph states with the same 
critical fidelity. 

Local Z-noise also corresponds to the thermal state of 
the Hamiltonian in Eqn. ([1]), and provided we set J; = 
A/2, the local error probability is the same at each site. 



1 e-^^ 



where the temperature T is encapsulated in the parame- 
ter (3 = 1/ (kBT). Given our proof of optimality, this can 
be phrased as a critical temperature. 



-A 



ln(V2 - 1) 



(4) 



which corresponds to one of the bounds found in 

as one expects since the argument of fl6'| also involves 



the breaking of cluster states into two-qubit states. As 
discussed in [13] , it should be possible to probe this tem- 
perature with existing experimental implementations. 




FIG. 6: (a) Take a general graph and bipartition it witli a 
single qubit on one side, (b) Apply local controlled-phase 
gates along all edges, (c) Apply unitaries y/K to change the 
central node of the GHZ state. Theorem[T]can now be applied 
to obtain upper bounds. 



Thresholds can also be derived for arbitrary graphs by 
first applying appropriate transformations. Select from a 
particular graph, G, a single qubit which has the min- 
imum number of nearest neighbors (minimum degree) 
Dmim which we take to constitute a partitioning of the 
graph. We now give a reversible procedure which is local 
to this partitioning to convert the state into a GHZ state 
which is LR across the partition. By applying controlled- 
phase gates along all edges that are not connected to our 
chosen qubit (remember that these gates commute with 
the Z-noise), we are left with a (Dmin + l)-qubit GHZ 
state as depicted in Fig. (HUb). However, the hnks across 
the partition are not the same as for an LR state. To 
account for this, we apply local unitaries y/TCi on our 
chosen qubit and one of its neighbors, performing the 
transformation depicted in Fig. [D^c). These local uni- 
taries transform the noise model from local Z-errors to 
an X-error on the original qubit, a K-error on the cen- 
ter of the GHZ state, and Z-errors everywhere else. The 
corresponding density matrix can be written in the graph 
state basis as 



p= Yl (l-pr^p'"^diag((l-p)^™+i-2-^p(l-p)^™-2-.,p^™+i-2-.,(i_p^^^ ^ (5) 

i6{o,i}°mi„-i 

where wj is the binary weight of j. These diagonal elements (j| p \j) are indexed by a number j G {0, 1}^""'^, where 
the J*'' bit indicates whether a Z-error has occurred on the z*'' qubit relative to the desired pure state \iPg)- The most 
entangled pair of qubits is given by j = 0, so we define 



diag ((1 - p)A„n + l^p(l_p)A„n^p I^n„n + l^ (l _ p)pI3„.„) 

(l_p)A„„+pl 



= r. — ^^D„„„ ^ • (6) 



r 



From this state, the entire state p in Eqn. (O can be re- constructed. Hence, if p2 becomes separable, purification 



7 



must be impossible, i.e. if 



P 



For -Dmin = 1, we recover the bound for the LR graphs as 
expected. However, LR graphs need not have -Dmin = 1; 
proving that this bound is not always tight. For the trian- 
gular and icosahedral graphs, these calculated threshold 
probabilities are p = 0.352 (which the GMPP appears to 
numerically saturate, thereby exceeding the DRPP) and 
p = 0.413 respectively. As Dmin ^ oo, p — > i. The 
bounds given by [l^ in this case also show that purifica- 
tion is impossible ii p > ^, but make no tighter claims. 



B. Maximally (Global) Depolarizing Noise 

In the previous subsection, we proved optimality of 
purification for local Z-noise. The local unitary equiva- 
lence of graphs can provide similar bounds for a range of 
other local noise models. It is now interesting to examine 
the case of correlated noise, and to derive bounds in this 
context. We will prove bounds for all graph states by 
initially restricting to linear cluster states. Choosing the 
noisy state to be purified as the maximally depolarized 
state of an 7V-qubit linear graph state \iPg), 
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we give an inductive proof which shows how to locally 
create pn{x) from pjq-i(x). If we take pN-i{x) and add 
an extra qubit to it, then the density matrix takes the 
form 

p'^_i(x) = diag(l + .T, 0, 1, 0, 1, 0, ... , 1, 0)/(2^-i + x). 

Upon application of a Z-rotation to the new qubit (the 
least significant bit of i), the zeros and non-zeros swap. 
Further Z-rotations on the other qubits can permute the 
position of the 1 + x term. Taking each of these with 
probability (1 — p)/2^~^, or p'pf_i{x) with probability p, 
we are left with the density matrix 

2^^diag(H-x,0,l,0,l,0,l...) 
+^diag(0,l,0,l,0,l...) 
This can be forced to take the form of Pn{v) by selecting 

2^-1 +x 



P 



and X — y. Consequently, the threshold value of a; = 2 
for the bipartite case holds for all TV, and the threshold 
fidelity is 



fi 



threshold 



2^ + 2' 



(7) 



By allowing an arbitrary two-qubit density matrix of the 
form p2 = diag(i, a, 6, ^ — a — 6), it is possible to prove 
that no better bound can be given by this method. 



In App.IbI we calculate the performance of the DRPP, 
which performs very poorly in the presence of correlated 
noise, being unable to purify if / < 1/3. For = 2,3, 
the GMPP manages to saturate this bound. However, for 
A > 3, we have been unable to find suitable repetitions 
of PI and P2 which purify if / < 1 /2^/2^ This is because 
the maximally mixed state of one set of errors (i.e. errors 
on just the 'A' qubits, with pure 'B' qubits), which has 
fidelity 1/2^/^, is a strongly attracting fixed point. In 
App. El we apply the results of [13] to show that for 
the closely related protocol of [2§| when the qubits are 
partitioned into two equally sized sets, / — 1/2^/^ is 
indeed a fixed point. 

The threshold fidelity for GHZ states is also given by 
Eqn. (I7|), since the above proof also holds for all states 
with a single edge across the bipartite division. The 
advantage is that the two potential fixed points in fi- 
delity (1/2 and l/2^~^) are different from those of a 
linear chain. Moreover, the smaller of these is below the 
threshold fidelity, and is trivially avoided. For example, 
for A^ = 5, purification of a; > 2.024 is possible using 
the GMPP. Given the anticipated asymptotic approach 
to the critical fidelity, the GMPP seems to saturate the 
bound of Eqn. 0. The threshold fidelity for GHZ states 
coincides precisely with that of [l^ , although our chosen 
parameterization of the state provides a more convenient 
condition for purification. 

It is now possible to follow an identical protocol to 
Fig.[6]in the case of maximally depolarizing noise to prove 
that all graphs are subject to the bound in Eqn. ([T]). This 
follows trivially from the observation that the application 
of controUed-phase gates and local unitaries \/Ki do not 
change the noise model, only the underlying graph. Once 
we have a GHZ state with a suitable partition, the above 
derivation applies. 



C. Local Depolarizing Noise 

In addition to local Z-noise, a range of other local noise 
models could be considered. One such model is where the 
type of local unitary that is applied is not known, but it 
occurs with probability p. This is equivalent to local 
depolarizing noise occurring with a probability 4p/3, 

Slip) = (i-p)p + |(A,pA, + y,pr, + z,pZ,) 



p 



4pi 
3 



^il,®Tr,(p). 



To derive a threshold, it is possible to follow a similar pro- 
cess to the previous subsection. Of critical importance to 
calculating a tight bound is the optimal selection of the 
two-qubit state to use across the partition. When con- 
sidering linear graphs, this state varies with the length 
of the chain. To illustrate this, we shall discuss the case 
of 3 qubits in more detail. The target density matrix, 
written in the graph state basis, takes the form 

p = diag(a, 6, c, b, b, b, b, b), 
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FIG. 7: Diamonds indicate the probability of a local depo- 
larizing error above which purification of a linear chain is 
impossible. The dashed line indicates a bound below which 
all such probabilities must lie. Stars indicate values for which 
purification can be achieved with the GMPP. 



where a, b and c are specified in terms of p. This can be 
divided into a probabihstic mixture of two components, 

p = diag(a,0,c,0,5,0,5,0) + diag(0,6,0,6,0,6,0,6) 

= diag(a, c, b, b) ® |0)(0| + ^diag(l, 1, 1, 1) ® |1)(1|. 

The second of these is a maximally mixed state on the two 
qubits that have an edge crossing the bipartite partition, 
perfectly connected to a third qubit (with a Z-rotation). 
It is always possible to prepare this, so our only condition 
relates to the creation of the first term, which is a two- 
qubit state diag(a, c, 6, b)/{a + 2b + c) perfectly connected 
to the third qubit. This is the form of two-qubit state 
that we choose to use, and purification is impossible if 
a/(a + 26 + c) < \. Upon evaluation of a, b and c, this 
gives 

27 - 126p + 156p^ - 6 V < 0, 

from which we find the threshold probability of 0.332. A 
similar process can be adopted for all linear graphs, and 
the results are presented in Fig. [T] In App. |B1 we de- 
rive the performance of the DRPP, which fails to achieve 
these bounds. For 2 qubits, the GMPP is the same 
as the protocol in j7|], and therefore achieves the two- 
qubit bound of 31.7% asymptotically. For 3 qubits, the 
GMPP comes close to matching these bounds, purifying 
at 33.1%. However, for A^ > 4, the strongly attracting 
fixed point of / = 1/2 ^^/^l occurs close to the threshold 
probability. This makes numerical analysis particularly 
challenging in these cases. 

These bounds exceed the bound for local Z-noise be- 
cause for the graph states, Z-noise forms a complete ba- 
sis, whereas the other errors do not. Consequently, when 
two errors coincide, they might cancel, and hence it may 
become slightly easier to purify the state. For example. 



(a) 
•- 

(b) 

•- 



FIG. 8: For the GMPP, purification is geometry independent. 
This means that the purification regimes for the two depicted 
states are identical. 



the purification condition for two qubits subject to lo- 
cal depolarizing noise is (1— p)^-|-p^/3 > | instead of 
(1 — p)^ > \ (.Z-noise), where the extra term comes from 
cancellation of coinciding errors. 

A general bound for a chain of arbitrary length can 
be obtained by demanding the conditions under which 
we can create a particular A^-qubit state, where only the 
first A^ — 1 qubits are noisy. Provided A^ > 3, we can 
create a chain of arbitrary length M > A^ by adding 
M ~ N qubits to the end of the chain, along with the 
required noise. This bound must be non-increasing with 
increasing A. Selecting A = 10, wc find that purifica- 
tion is impossible for all chains of ten or more qubits if 
p > 0.347. Consequently, the threshold probability must 
tend towards a constant, as observed numerically for the 
GMPP in [l3|. The optimal protocol must also tend to 
a constant because of the constant lower-bound provided 
by the DRPP (App. E]). 

In comparison to linear or GHZ states, the case of 
the Steane [[7,1,3]] code is slightly more involved be- 
cause we need three two-qubit states to cross the bipartite 
split. However, all three become separable at the same 
threshold probability of p = 0.403. Numerically, we find 
that the GMPP becomes trapped in the fixed point of 
/=l/24. 



D. Analysis of the GMPP 

To date, evaluation of the performance of the GMPP 
has resisted analytic techniques, and has instead relied 
on numerical evaluation, as we have used in this section. 
While we have observed a close relationship between the 
GMPP and the analytic bipartite purification protocol of 
[29I [30j , the direct proof of any connection still remains 
an open problem. However, we are able to make progress 
in proving the purification regime for the GMPP in cer- 
tain special cases. Given this, it becomes interesting to 
compare the performance of the two purification proto- 
cols, the DRPP and the GMPP. Earlier in this section, 
we have provided several examples where the GMPP out- 
performs the DRPP. In this subsection we will present 
our analysis of the GMPP, and construct an example 
in which the GMPP is provably sub-optimal, being out- 
performed by the DRPP. This yields an interpretation 
as to why the GMPP can get trapped in fixed points for 
certain error models. 

We proceed by realizing that the sub-protocols of the 
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FIG. 9: Consider purification of the graph in (a) where gray 
circles denote qubits with ^-errors with probability p, and 
black circles are pure qubits. For the GMPP, purification 
of parts (a), (b) and (c) are identical, (b) provides an upper 
bound of p = 30% for the GMPP, and (c) proves that this can 
be achieved. The DRPP can always purify (a) for arbitrary 
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FIG. 10: Local operations can convert local Z-noise on graph 
(a), which is LR, into a two-colorable graph (b). The GMPP 
can purify this if it can purify (c). However, this is not the 
same constraint as the ability to purify a single copy of the 
three-qubit state. 



GMPP depend only on the diagonal elements of the den- 
sity matrix, and not the underlying geometry, other than 
the numbers of qubits of each color in the two-colorable 
graph. This means that the purification regimes of the 
graphs depicted in Fig.[8]are identical. Hence, one can re- 
strict to considering purification of the graph in Fig.[Sl^b). 
Further, if we assume that the noise is not correlated be- 
tween the pairs, and that it is identical for each pair, 
then the GMPP is exactly the same as the DRPP, ex- 
cept that purification of the pairs occurs in parallel in- 
stead of independently, leading to the observed reduction 
in purification rate. Given that we can derive the perfor- 
mance of the DRPP, we can deduce the performance of 
the GMPP in this case and, consequently, in the case of 
more complex graphs such as Fig. [Hfa). One such exam- 
ple of noise is local Z-noise, instantly proving that the 
GMPP is optimal for local Z-noise on LR graphs. 

The geometry independence of the GMPP means that 
it can purify all two-colorable graphs with local Z-noise 
with the same critical probability. This includes graphs 
such as Fig. [ijd), for which our analysis gives a thresh- 
old probability oi p = 0.352. We interpret the failure 
to saturate this bound as the geometry independence of 
the GMPP causing it to become trapped by local fixed 
points. 

In Fig. [21 we consider purification of a linear graph 
where two of the qubits have a Z-error with probability 
p, and the other two qubits are pure. Given the geometry 
independence of the GMPP, its critical probability is the 
same as for Fig. IH^b). As we have seen in Sec. IIV Al 
purification of this graph must be impossible if p > 1 — 
1/V2. Consequently, the GMPP cannot purify Fig.[9l^a) 
if p > 1 — 1/ \/2- A similar manipulation to Fig. [8] yields 
Fig.lHc), which can also be used to show that purification 
below this critical probability is possible (application of 
purification to the pure pair leaves it pure, so we only 
have to purify the noisy pair, for which the GMPP is the 
same as the optimal two-qubit protocol). 

Now consider applying the DRPP to the original chain 
(Fig. inia)). Each of the three edges to be purified can 
always be purified - one is already pure and the other 
two only have two diagonal elements to their density ma- 



trices, so can always be purified. Hence, the state in 
Fig. IHa) can always be purified. This proves that the 
GMPP is sub-optimal and that in some circumstances 
the DRPP out-performs it. 

There are certain pitfalls associated with this analysis 
of the GMPP that we will illustrate with an example. 
Consider the purification of the triangular graph with lo- 
cal Z-noise, for which our upper bound has predicted the 
impossibility of purification if p > 0.352. The state can 
be transformed into a linear graph by local operations 
where the errors are now Z,Y, Z , and since this graph is 
two-colorable, the GMPP can be applied. Numerically, it 
appears to saturate our bound. We can now consider the 
purification of the LR state depicted in Fig. [TUTa). With 
local Z-noise, we already know that purification is impos- 
sible if p > 0.3. However, we can apply local operations 
to form a two-colorable state, and subsequently apply the 
geometry independence of the GMPP to see that purifi- 
cation of this state is equal to the parallel purification of 
two of the triangular graphs, and hence it appears that it 
should have the same purification regime. Clearly, there 
is a discrepancy. This is resolved by observing that the 
performance of the purification in parallel is not identi- 
cal to two independent purifications in this case because 
there is an asymmetry between the PI and P2 protocols 
i.e. the parallel application requires PI' = PI (g) PI and 
P2' = P2 (g) P2, where PI and P2 are the protocols on 
the triagle, whereas the GMPP applies PI' = PI (g) P2 
and P2' = P2(g)Pl. 



V. PURIFICATION OF VALENCE BOND 
STATES 

Using a combination of the DRPP and Theorem [U we 
have proved that for local Z-noisc, and some other types 
of local noise, there is a threshold error probability below 
which purification is possible, and above which purifica- 
tion of graph states is impossible. However, there is no 
need to restrict to graph states. By making use of the va- 
lence bond formalism, we will now construct an example 
of a state which is not a graph state, but can be optimally 
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purified by the DRPP. 

A general valence bond state [H, HH, [H, [111 can be 
described by employing a D-dimensional maximally en- 
tangled state between each nearest-neighbor of a graph. 
Each local party then projects down to a d-dimensional 
system with a specific projector. If we allow D ^ 2^/2, 
then any 7V-qubit state can be described by this formal- 
ism [311 . The class of translationally invariant states 
can be described efficiently, using a fixed D. These D- 
dimensional maximally entangled states can be formed 
from log2(£') Bell states. We are solely interested in 
constructing a simple example to demonstrate the gen- 
eral properties. As such, we shall restrict to D = d = 2 
and to a linear graph of 3 qubits. This contains all the 
essential properties of valence bond states and, conse- 
quently, we expect that generalizations will follow in a 
straightforward manner. 

The purification protocol follows the concept of the 
DRPP, as already outlined. Our initial state is described 
by two maximally entangled states joined by a single 
projector Pq, acting on qubits 2 and 3. 

IV'initial) = 1 ® Po ® 1 !<?!>) 12 10)34 

If local Z-noise affects our state, then we apply Z- 
measurements to all the qubits apart from a single pair, 
which should retain some entanglement. Since the noise 
commutes with the measurement, it suffices to describe 
what happens to the pure state, 

|V2)|0> - il®Po®|0>(0| fEl*)l*)j 

= ^l(SPo®t(^\i)\i)^ |0>|0), 

where we have assumed outcome |0)(0| from the Z mea- 
surement. Using many copies, the state \ip2) with local 
Z-noise can be purified to a two-qubit maximally entan- 
gled state, 10). We repeat this for each edge of the graph, 
and the pure state |'0initiai) is recovered by applying local 
projectors Pq to each vertex. We can express the projec- 
tor Pq as 

Pa = Y^ aJ,^k \j) {i\ {k\ 

so that a'' constitutes a. d x d matrix, {j \ \i) = aj,ik- 
Note that this definition does not coincide with the stan- 
dard matrix product state definition of these matrices 
(where (j| \i) = ak,ji). 

For the optimality proof, we need to start with the 
state 1^2) and show how to reconstruct IV'initiai)- We do 
this by locally introducing a maximally entangled state 
10), and applying a projector Pi. 

l(g)Pl(g)l 1^2) 10) = llZ-initial) • (8) 
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FIG. 11: The valence bond state that we wish to purify is 
generated from two maximally entangled states and a 
projector P applied between them. 

If Z- noise is present on \'ijj2), then it must reappear on 
iV'initiai) when we apply Pi, 

(l®Pl®l)-(l®Z0l0l) IV'2) 10) = a®Z®l|V'initial)- 

(9) 

Pi can be described analogously to Pq, 

which, through Eqn. ([5]), allows us to show that 

Subsequent expansion of Eqn. ^ enables the derivation 
of a simple condition for when the reconstruction can be 
performed, and hence when the optimality proof holds, 

[f3\Z]a°^0. (10) 

This does not hold for all valence bond states, but we 
can construct examples when it does. In the case where 
a'^ and are invertible, we find that a'^(a^)~^ must be 
diagonal. Applying the optimality proof on both edges 
of the graph provides a symmetry between the elements 
aj^ik and aj^ki- This leads to a final form of the projector 

— ^j^-jTj^ — ai,oi ai.io / 

In the special case of ai^n — e'^^i"'"^^)^ ai^io = e*^^, 
0^1,01 = e'^^ and ao,ii = ao.io = ao,oi — 1, we recover 
the weighted graph states [H, |3l] and the cluster state 
(6*1 = O2 = it). Weighted graph states are identical to 
the graph states that have been discussed so far except 
that to construct the state, instead of a controlled-Z gate 
between nearest-neighbors initially in |-|-), a controUed- 
phase gate of arbitrary phase is used. This means that all 
the relevant actions continue to commute with Z-errors, 
and we recover trivially the previous optimality proof, 
providing a useful verification of these results. Since the 
weighted graph states have an exponentially decreasing 
localizable entanglement length i35i] , the general so- 
lutions, as described by Pq, are expected to have finite 
localizable entanglement length. 

VI. CONCLUSIONS 

In this paper, we have described a method which 
proves that certain noisy multipartite states cannot be 
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purified. For the case of LR states subject to Z-noise 
with probabihty p, we have been able to show that for 
p > 30%, purification is impossible, and that all other 
states can be purified i.e. we have demonstrated opti- 
mality of the purification protocol. 

Numerical evidence indicates that the GMPP [Til . [T3 | 
is optimal (in terms of the states that can be purified) 
for a large range of different errors, including local er- 
rors such as Z-noise (for which we have proven optimal- 
ity), and non-local errors such as maximally depolarizing 
noise. There are cases where the GMPP is not optimal 
and in all such cases, we have observed that the protocol 
gets trapped by strong local attractors with fidelities ^ 
and where the two-colorable state has n qubits of 
one color, and m qubits of the other. These can often be 
interpreted as being due to the geometry independence 
of the GMPP. This also coincides with the results that 
can be derived for the bipartite purification protocol (ap- 
plied to systems of arbitrary dimension) in [2^, [s^l ■ The 
purification of the noisy graph in Fig.[HJa) provides proof 
that this is a real phenomenon, and not just an artefact 
of finite computational resources - the GMPP has a crit- 
ical probability for purification, but other protocols can 
always purify the graph. 

Not only have we demonstrated the optimality of pu- 
rification of some graph states, but have also provided 
an example of a valence bond state which can be opti- 
mally purified, thus demonstrating the general utility of 
our method. 

In forthcoming work, we shall examine how our upper- 
bound method can be applied to the situation where the 
gates used during the purification are also faulty, which 
is a major advantage of our constructive approach. This 
has potential implications for upper bounds on fault- 
tolerant thresholds. Interesting extensions of this work 
could involve taking what we have learnt about the per- 
formance of the GMPP and trying to improve it. In 
particular, we have demonstrated that one should take 
into account both the geometry of the state and the noise 
model when constructing a purification protocol, not just 
the noise model (in the case of the GMPP) or the geom- 
etry (in the case of the DRPP). A combined approach, 
based on the stabilizers of the state, but allowing asym- 
metry between the terms, appears to be the most sensible 
approach. 
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APPENDIX A: ANALYSIS OF LR CLASS 

Given our rather abstract definition of the LR class, it 
is a worthwhile exercise to verify that it is a non-trivial 
class. 

1. Graphs of Maximum Degree 3 

Let us consider all possible graphs which have a maxi- 
mum number of nearest neighbors (degree) equal to 3 or 
less. If one of the qubits in the graph has degree 1, then 
there exists a trivial partition that shows that the state 
is LR. 

Arbitrary graph 

i ' 

/•■-; 

Partition 

All other graphs must contain a loop. If we take the 
smallest loop in the graph, then the qubits in this loop 
already have at least two neighbors. We start by tak- 
ing the triangle, a loop of 3 qubits. With no additional 
connections, this state is non-LR. 

V 

If we add one extra qubit, then by connecting it to two 
or three qubits in the triangle, the state is non-LR. 

When the extra qubit is only connected to two qubits, 
there are two further links that could be added to any 
arbitrary structure. If they are not connected to the same 
qubits, these two links provide an LR partition. 

The only remaining structure is where all three links from 
the triangle are connected to some arbitrary graph, but 
are not incident on the same qubit. These extra links 
must provide an LR partition. 

We now continue this argument to loops of 4 qubits. 
Any structure which we add that generates a triangle 
has, of course, already been dealt with. This only leaves 
four examples which are non-LR. 

Finally, for larger loops, there is no way of creating a 
graph that is non-LR without forming smaller loops first. 
Consequently, for graphs of maximum degree 3, there are 
only 8 graphs which are non-LR. It can be verified that 
the latter are locally equivalent to LR states. 



12 



N 


LR states 


All graph states 


3 


7 


8 


4 


53 


64 


5 


788 


1024 


6 


22204 


32768 


7 


1148781 


2097152 



TABLE I: Comparison of the number of LR states of A'" qubits 
and the total number of graph states. Note that we have 
included the completely separable state in the set of LR states. 



2. Size of LR Class 

It is important to answer the question of how large the 
set of LR states is. The set of all graph states of qubits 
consists of 2(") elements, including all possible isomor- 
phisms and graphs which can be separated into two or 
more unconnected sub-graphs. We can easily generate 
LR states from the (A^ — l)-qubit graphs by adding an 
extra qubit and connecting it to any single qubit from 
the previous graph. This state must be LR because the 
qubit that we have just added constitutes one such par- 
titioning. There are at least (A^ — 1) such states. 
Despite being a small fraction (~ of all graph 
states, they certainly form a significant class in their own 
right. Moreover, there are further examples to be added 
into the class of LR graph states, but accurate enumera- 
tion is a combinatorial challenge. In Table [T] we present 
the number of LR states for graphs of up to 7 qubits. 

Now that we have provided a lower bound, is it possible 
to either tighten this bound, or give an upper bound? 
Let us consider all possible graph states of A^ qubits. 
We shall select a specific partition q < N — q qubits. 
The probability, Pq, that local reconstruction is possible 
across this boundary is given by 



Pq = 



b=0 



)(V)^! 



2<}{N-q) 



This is a result of requiring that there are q{N — q) pos- 
sible bonds across the partition, which could either be 
bonded or not. Of these 2''^^"'^ combinations, only those 
with no more than a single bond from each qubit make 
the state locally reconstructible. For b bonds across the 
partition, we have to choose them from q on one side 
and A^ — q on the other side. Finally, the ordering of the 
choice on one side of the partition is important, hence 
the 6!. This can be expressed in terms of the confluent 
hypergeomctric function. 



Pq 



= 2«('J-^)(-l + {-iyU{-q, l-2q + N,-l)). 



There are (^) ways that we could have chosen a parti- 
tion of q qubits. Each of these has the same probability 
of giving local reconstructibility, but we must make sure 
we do not over count the cases where there is more than 



one partitioning for the same graph state, 

n=l \ / 

Similarly, we need not be restricted to a specific but 
must avoid over-counting, 

LJV/2J 

To simplify this expression, we can take the smallest 
(largest) value of Pq and assume that all Pq have this 
value, thereby lower (upper) bounding P. 



P > 



LW/2J 

E 

n=l 



N/2 



-'mill ~ -I minj 



N/2 



The bounds will occur for q ~ I and q = N/2. In the case 
of q — 1, Pq — {N — l)/2^^^ and hence the fraction of 
LR states is estimated to be N^ /2^ . To see that this is 
an upper bound, we write that Pq ~ 29('J-^)0(A^9) and 
P ~ Pq (^) Y ' having assumed that a particular value 
of q is going to give the required bound. The ratio for 
successive values of q is therefore given by 



PI, 



N -q 



\q-l 



■ 0{N) ■ 229- 



-N-l 



Hence, each successive value of P must be smaller given 
the overpowering nature of the exponential 2~^. There- 
fore the bound which we have just derived is an upper 
bound. Combining this with our existing lower bound, 



A^-1 ^ \LRn\ 



\G 



< 



N 



N3 



APPENDIX B: THE PERFORMANCE OF THE 
DIVIDE AND REBUILD PROTOCOL 

In the body of the paper, we focused on calculating 
upper bounds for certain types of noise, and comparing 
them to a numerical analysis of the performance of the 
GMPP. We have also analyzed the DRPP in the case of 
local-Z noise, when we can show that it is optimal. In 
general, we do not expect this protocol to be optimal, but 
it is still useful because we can analyze its performance, 
and use it to place a lower bound on the performance of 
any optimal protocol. 



1. Maximally Depolarizing Noise 

We can readily show that the DRPP does not adapt 
well in the presence of some correlated noise (of course, X 
and y-errors can be represented as correlated Z-errors, so 
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we have optimality in some special cases). For example, 
we can take the case of the maximally depolarized state 
of N qubits, 



IN J_rr 



Pn = 



and consider \iPg) to be the A^-qubit linear cluster state. 
When we perform measurements on this state, we reduce 
it from pn{x) —t pN^i{x/2) p2{x/2^^^). The bipar- 
tite state can be purified if 



1 



2"-3 1 
X ^ 1 



and hence the fidelity goes as 



3 3-2^-1 

which tends to a fixed value of 1/3, whereas genuine mul- 
tipartite purification protocols can purify states exponen- 
tially decreasing fidelities [l^ . 



2. Local Depolarizing Noise 



to the final qubit. 



n=0 
DG-l~2n 



2n 



m=0 ^ 

i(l + (l-V3)^--i) 



If g{n) is the probability that n errors occur on the Bell 
pair due to its local errors (i.e. neglecting the connected 
qubits), then after measurements there is no error with 
a probability 

g{Q)q^+g{l)q{l-q)+g{2){l~qf. 

No errors occur on the Bell pair if no errors truly oc- 
curred, (1— p)^, or pairs of errors cancel {YY , ZX, XZ), 
p^/3. Hence, g{0) = (1 — p)^ +p'^/3. Similarly, we find 
that 5(2) = ^g{l) = 1 - g{0)/3. Substitution of these 
yields the polynomial 

a;2Do ^ 2X^« + 1 - 1 - 

where x = 1 — 4p/3 and p is the critical error rate below 
which this protocol will perfectly purify the state. In the 
case of the Steane code, Dq = 3, and hence Pcrit = 

0.16. 

As stated in the body of the paper, existing multipartite 
purification protocols improve upon this probability. 



We would also like to demonstrate the performance of 
the DRPP when we do not know what the type of lo- 
cal error is. We proceed by assuming that the noise is 
the most destructive type of local noise considered here. 
This is simply Z-noise, because these errors form a com- 
plete basis for the state, whereas other types of error 
need not. Errors occur with probability p, with an equal 
likelihood of them being either X, Y or Z. Making a 
Z-measurement hence causes the propagation of an error 
to another qubit in 2/3 of the cases. The simplicity of 
the protocol, however, continues to help us, enabling the 
calculation of the critical error probability, Pcrit assuming 
that the graph has a maximum degree of Dq- Note that 
it is necessary to assume that the graph is two-colorable, 
which means that an error can only propagate to one of 
the two qubits in the Bell pair to be purified. 

Consider a single Bell-pair, which is the one which we 
aim to measure towards and purify. Each qubit is at- 
tached to Dg — 1 other qubits in the graph. We will be 
able to purify the state if the probability of error on the 
Bell pair, after measurement of the other qubits, is better 
than i . The probability of an error occurring on a single 
qubit is p, and it is equally likely to be X, Y or Z. Let q 
be the probability that all the qubits attached to one of 
the pair give no errors on the qubit they are attached to 
after Z-measurements. This is caused by no errors occur- 
ring, Z errors occurring (which do not get transmitted) 
or pairs of X or y errors, which cancel when transmitted 



3. Fully Connected Graph States 

As an example of the application of the DRPP, we will 
examine the purification of fully-connected graph states. 
Note that these graphs are not LR, although they are 
locally equivalent to GHZ states, and hence have been 
indirectly included in previous discussions. We shall start 
by considering the triangular configuration of Fig. [5jc) . 

Given the small size of the triangle, examination of 
the different local error combinations is tractable. There 
are 3"^ = 27 different combinations of local errors. Of 
these, 12 can be turned into independent local .Z-noise on 
the 3-qubit chain using the equivalence of graphs under 
local unitaries. These are the combinations XXY, XYZ 
and ZZY and their permutations. Two further cases of 
particular interest are local X-noise and local F-noise. 
Our protocol of measuring a single qubit (in the X-basis 
in this case) can tolerate an error probability of p < i for 
X-noise. This is because pairs of errors obey identities 
such as ^2X3 = Xi, and hence an A"-measurement on 
qubit 1 commutes with these errors. 

The ultimate realization of this is in the case of in- 
dependent local y-noise, since a y-error is the same as 
correlated Z-errors on all 3 qubits. Hence, all the Y- 
errors are the same, and a single K-measurement would 
eliminate all of them. However, a single y-measurement 
leaves all the other qubits in a separable state - all the 
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bonds are broken, not just those connected to the qubit 
that is measured. Instead, we must perform some other 
measurement. A Z-measurement, for example, allows 
purification for all errors since one of the two diagonal 
elements is always larger than i. 

The results derived for a triangular graph for both local 
X-noise and local F-noise hold for all fully connected 
graphs. 



qubit, and Z-rotations on all the other qubits. Since the 
independent X errors commute with this measurement 
process, this will eventually reduce to the triangle, which 
we have already solved, and hence p < ^ is the criterion 
for purification. Note, however, that when reducing to 
a two-qubit state, our density matrix has 3 diagonal ele- 
ments, not 2 as in the case of Y errors, so purification at 
p > i is not possible with this protocol. 



\Gn)^ 



1 



where Wi is the binary weight of the number i. In the 
case of F-noise, a single y-error is the same as correlated 
Z-errors across all qubits, so all K-errors are the same, 
and the noisy state only has two diagonal elements, the 
first of which is the probability that an even number of 
errors occurred, 



L7V/2J 
ri=0 



= i ((1 - 2P)^ + 1) 

Upon performing a Z-measurement on one of the qubits, 
we are left with the state \Gm-i) and the diagonal ele- 
ments must be the same as they were before. This can 
be verified by writing 



\Gn) 



1 



E 

3=0 



(io) + (-ini))(-i)("^)|j) 



(Bl) 

Note that (-1)'"^ \j) = Z®^ \j) if j is an iV-bit binary 
string, which means that if we get the |1)(1| result, the 
corrective unitary is a Z-rotation on every qubit. This 
reduction to fully connected graphs continues to the case 
of IG2), which we know can always be purified if it only 
has two diagonal elements since one element is always 
greater than i. 

In the case of X-noise, we must verify that an X- 
measurement on a single qubit reduces \Gn) to \Gn~i)- 
After applying to the state in Eqn. (|B1|) . we are 

left with a sum over binary strings of even weight. To 
convert this into IGat-i), we apply a Hadamard on one 



APPENDIX C: FIXED POINTS OF BIPARTITE 
PURIFICATION 

In [2^, a bipartite purification protocol was proposed 
for which subsequent analysis [30| showed that its behav- 
ior is very similar to that of the GMPP, with the added 
advantage that we can calculate its performance. As an 
example, let us consider maximally depolarizing noise ap- 
plied to a linear state of N qubits. For simplicity, and 
consistency with [30[, we shall define D = 2^/^ assum- 
ing that there are N /2 qubits either side of the bipartite 
split. The diagonal elements of the initial density matrix, 

\^^^ indicate the errors k and j either side of the parti- 
tion. After n iterations of the protocol, the mmormalized 
outcome is 



\j - 



D-1 

E' 

k'=0 



D-1 

E' 

Lk"=Q 



-, (2") 



For maximally depolarizing noise. 



A 



(0) 
kj 



1 



+ SjSk 



Following this through and rcnormalizing, wc find that 

(n) _ jx + Py'" + {D~ 1)X^" 



A, 



00 



D[{x + DY-^ -h (£)- 1)D2"]' 



For X — D, Xqq is independent of n, and hence we have 
found a fixed point of the protocol. The implication is 
that there is a critical fidelity of 1/2^/^, below which the 
maximum achievable fidelity is 1/2^/^, and above which 
complete purification is possible. 
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